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ABSTRACT 

A  study  of  dissociation  reaction  Of  a  diatomic  gas  has 
been  performed  for  the  pressure  history,  a  pressure  rise 
and  subsequent  decay,  common  to  hypersonic  flow  about  a 
body.  For  this  Simplified  case,  several  interesting  phe¬ 
nomena  which  occur  have  been  explored.  A  characteristic 
reaction  time  has  been  identified;  the  areas  in  which  dis¬ 
sociation  or  recombination  is  the  faster  reaction  have  been 
specified;  the  areas  in  which  equilibrium  isentropic  expan¬ 
sion  leads  to  dissociation  or  recombination  are  found;  the 
relative  importance  of  the  irreversibility  associated  with 
nonequilibrium  relaxation  has  been  defined. 
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NOMENCLATURE 

A  Dissociation  reaction  rate  Constant 

B  Recombination  reaction  rate  constant 

Cv  Specific  heat  at  constant  volume 

Dissociation  reaction  rate 

k  Recombination  reaction  rate 

r 

M  Molecular  weight 

N  Reacting  species 

p  Pressure 

q  Velocity 

R  Gas  constant 

t  Time 

T  Temperature 

X  Nonreacting  species 

a  Mass  fraction 

0  Characteristic  temperature 

t  Characteristic  reaction  time 

«  Temperature  exponent  in  reaction  rate  equation 

Subscript 

eq.  Equilibrium 

f  Final 

i  Initial 

oo  Undisturbed  free  stream 

Superscript 

—  Average 

(N)  Mole  concentration 
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1.  INTRODUCTION 

The  flow  about  a  body  moving  at  supersonic  speed  is  characterized 
by  a  pressure  rise  through  a  bow  shock  wave  and  an  expansion  about  and 
behind  the  body  to  ambient  pressure.  The  details  of  the  pressure  pattern 
depend  on  the  exact  body  shape  but  the  basic  compression  and  expansion 
features  are  always  present. 

At  hypersonic  speeds,  where  the  kinetic  energy  of  the  flow  relative 
to  the  body  is  Of  the  Order  of  the  dissociation  energy  of  the  air  molecules, 
such  reactions  can  be  expected  to  take  place.  For  a  pure  diatomic  gas, 
such  as  nitrogen,  only  the  dissociation  reaction  will  take  place  while  air 
will  undergo  several  reactions  involving  both  oxygen  and  nitrogen.  All 
these  reactions  proceed  at  a  finite  rate  so  that,  depending  on  the  size  of 
the  body,  the  ambient  density,  and  other  characteristic  parameters,  the 
chemical  reactions  may  stay  essentially  in  equilibrium  with  the  pressure 
changes  or  may  badly  lag  behind  the  pressure  changes.  In  some  cases 
the  chemical  reaction  are  so  slow  that  they  may  be  considered  frozen  to 
the  Scale  of  the  body,  but  if  the  condition  continues  long  enough  they  will 
eventually  approach  equilibrium.  A  nonequilibrium  condition  may  only  be 
considered  frozen  with  respect  to  some  particular  scale.  In  many  cases, 
flow  field  parameters,  such  as  velocity  and  pressure,  can  be  considerably 
effected  by  nonequilibrium  chemical  reactions  so  that  both  the  flow  field 
and  chemical  equations  must  be  solved  simultaneously. 

In  this  report,  a  diatomic  gas  subject  to  a  pressure  rise  and  decay 
separated  by  a  region  of  uniform  pressure  is  considered  to  bring  out,  at 
least  qualitatively,  the  interesting  phenomena.  The  rate  of  the  chemical 
reactions  taking  place  are  considered  in  order  to  specify  the  character¬ 
istic  times  and  the  parameters  that  control  these  times.  The  general  flow 
field  is  examined  to  discover  other  interesting  phenomena  which  are  then 
considered  in  some  detail.  The  pressure  rise  and  decay  are  both  con” 
sidered  to  be  discontinuous  step  functions  and  the  chemical  reactions  are 
considered  to  occur  under  conditions  of  constant  pressure.  The  time  for 
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chemical  reactions  to  occur  under  these  idealized  conditions  can  be  used 
to  determine  the  relative  speed  of  the  reactions  compared  with  the  actual 
rate  Of  change  of  the  pressure. 


2.  REACTION  RATE  ANALYSIS 
The  dissociation  reaction  for  a  diatomic  gas  can  be  written 


N2  +  Xt=;N  +  N  +  X 


(1) 


The  rate  equation  for  this  reaction  is 


^  =  2  kd(N2)(X)  -  2  kr(N)2  (X) 


(2) 


If  a  is  the  mass  fraction  of  N  then 


(N)  = 


(N,)  -  (X)  =4, 


*h' 


N, 


(3a) 


where 


l  a  .  1  -  a 


(3  b) 


being  an  average  molecular  Weight.  Using  the  perfect  gas  relation 


jEL  =  _E_ 
M  RT 


(4a) 


where  R  is  the  universal  gas  constant.  Applying  relation 


N 


1  +  Q 


(4b) 


since  =  2  Mj^.  Substituting  relations  (3)  and  (4)  into  (2)  the  rate  relations 
can  be  written  in  terms  of  mole  fractions,  pressure,  and  temperature. 

2  ,  _  n2 


dt  =  (l  "  q)  R?T  kd  ’  1  +  a  (l£r)  kr 
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(5) 


For  many  diatomic  gases 


kd  =  AT‘“  i‘0/T,  kr  =  BT"W  (6) 

For  Oxygen  and  nitrogen  there  is  fairly  general  agreement  on  the  values  of 
the  constants  in  these  equations.  Values  of  these  constants  which  seem 
best  based  on  a  current  survey  of  the  literature  are 


B 


OJ 


Oxygen  1( 1O20) 


cm 


mole  sec°K  u 


5(1018) 


cm 


,2 _ o*>~<a) 

mole  sec  K 


1.5  58,600°K 


Nitrogen  6( 1O20) 


3{1019) 


1.5  113,260  K 


These  are  the  numerical  constants  which  will  be  used  where  required  in 
this  report.  If  these  expressions  for  kd  and  kr  are  substituted  into  the 
rate  equation 


da 

dt 


ts  m  *  1 


RT 


+w 


,-0/T  _  4  a2  BP2 

1  1  +  a  r2t2+w 


(7) 


Equilibrium  conditions  occur  when  da/dt  =  0 


2 

or 


A  RT  ,-0/T 
41  “P-  1 


(8) 


or 

a2  A  1  -0/T 

1  -  a  '45  p  i 

Equation  (8)  is  consistent  with  Lighthill's  equilibrium  equation  for  an  ideal 
dissociating  gas  (1)  and  Equation  (7)  with  Freeman's  dissociation  rate 
law  (2)  In  these  references  the  term  A/4B  is  called  a  characteristic 
density. 
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The  next  step  in  considering  Equation  (7)  will  be  to  attempt  to  find 
a  characteristic  time  of  reaction  so  that  the  speed  of  reaction  may  be  easily 
estimated  It  is  also  useful  to  consider  which  parameters  effect  the  rela¬ 
tive  speed  of  the  dissociation  and  recombination  reactions.  A  simple 
answer  to  this  question  based  on  the  fact  that  the  dissociation  reaction, 
the  process  of  the  molecules  breaking  up  into  atoms,  requires  a  two-body 
collision  and  the  recombination  reaction  requires  a  three-body  collision 
is  not  sufficient.  While  the  two  molecule  collision  may  seem  more  prob¬ 
able,  this  does  not  mean  that  the  dissociation  rate  is  the  more  rapid.  At 
equilibrium  dissociation  and  recombination  are  proceeding  in  equal  rate, 

The  effective  dissociation  or  recombination  rate  is  the  rate  at  which  the 
concentrations  proceed  towards  equilibrium  if  the  initial  condition  is  not 
an  equilibrium  one.  This  rate  is  the  difference  between  the  dissociation 
and  recombination  rate  which  is  zero  at  equilibrium.  A  first  app r oxim ation 
to  both  the  effective  dissociation  and  r ec ombination  rate  can  be  obtained  by 
considering  the  rate  of  change  of  the  effective  rate,  with  species  concen¬ 
tration  for  an  adiabatic  change  at  constant  pressure,  evaluated  at  equili¬ 
brium  concentration.  In  this  approximation,  the  effective  recombination 
and  dissociation  rates  are  the  same  at  equal  changes  of  species  concen¬ 
trations  from  the  equilibrium  but  in  opposite  directions.  In  order  to 
distinguish  which  of  these  is  the  fastest,  the  next  higher  derivative  of  the 
rate  of  change  of  effective  rate  with  concentration  must  be  considered. 


The  first  approximation  to  the  rate  may  be  written 


teristic  time  t  can  now  be  defined  as 


and  Equation  (9)  can  be  integrated  to  give 
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(10) 


T 


a  in 


a*  =  a 

£..  ,  e,q 


a  -  a 


eq 


Where  a  is  the  initial  composition  and  the  final  composition. 


This  relation  shows  that  near  equilibrium,  where  this  first  order 
logarithmic  term  predominates,  all  reactions  should  follow  a  straight  line 
on  a  semi-log  plot  using  the  coordinates  of  Equation  (10)  To  determine 
how  well  the  actual  reaction  follows  this  line  at  large  deviations  from  the 
equilibrium  concentrations,  Equation  (7)  will  be  integrated  for  different 
initial  conditions. 


r  can  be  related  to  the  other  variables  by  differentiating  Equation  (7) 
and  evaluating  at  equilibrium  conditions. 


(11) 


If  equation  (11)  for  1/t  is  used  to  nondim ens ionaliz e  Equation  (7)  the  result 

is 


This  expression  can  be  integrated  if  T  is  known  as  a  function  of  a. 
The  relation 


can  be  obtained  by  expression  enthalpy  as  a  function  of  temperature  and 
composition  and  then  solving  for  the  temperature.  For  the  case  of 
inviscid  constant  pressure  reaction  h  is  a  constant  along  a  stream  line 
and  Equation  (13)  may  be  differentiated  to  give 


i  dT 
‘  T  Jq 


R 

C 

h  aR  + 


5“0  ‘  C~ 
v  v 


R 

C- 


1  +  (1  +  a) 


(14) 


Equations  (13)  and  (14)  may  be  used  in  the  integration  of  Equation  (12)  to 
give 


The  indicated  integration  of  Equation  (15)  has  been  carried  out  numerically 
for  R/Cv  =  1/3.,  cj  =  1.5,  AR0/Bp  =  10®,  10®,  10*®,  and  h/C^G  between 
0.1  and  0.5.  R/Cv  ~  1/3  is  the  value  suggested  by  Lighthill  and  probably 
the  best  compromise  if  the  variation  of  vibrational  specific  heat  with 
temperature  is  not  considered.  The  quantity  AR0/B  has  units  of  pressure 

7 

and  varies  with  the  gas  considered.  For  oxygen  it  is  2.4  (10  )  atm  and  for 

7 

nitrogen  it  is  4.7  (10  )  atm  so  the  chosen  values  are  in  the  range  of 
practical  interest  for  these  gases.  The  results  of  this  integration  are 
shown  in  Figure  la,  b,  c  plotted  to  the  coordinates  suggested  by 
Equation  (10). 

The  purpose  of  this  plot  is  to  compare  the  prediction  of  the  simple 
logarithmic  Equation  (10)  with  the  exact  solution  given  by  Equation  (15) 

The  departure  of  the  curves.  Equation  (15),  for  the  various  values  of  the 
parameters  from  the  straight  line,  Equation  (10),  shows  the  difference 
between  these  two  relations.  Before  comparing  in  more  detail  an 
explanation  of  the  integration  of  Equation  (15)  is  needed.  Since  the  time 
to  reach  equilibrium  is  infinite,  as  shown  by  either  Equation  (10)  or  (15), 
the  time  to  reach  some  final  composition,  near  equilibrium,  must  be 


-6- 


considered.  For  this  calculation,  a.  is  taken  as  far  from  equilibrium  as 

possible;  =  0  for  all  dissociating  cases  and  <l  =  1  for  all  recombination 

cases  except  those  cases  for  which  h/C^O  is  too  low  to  provide  enough 

energy  for  complete  dissociation.  For  these  cases  cl  is  taken  as  that 

value  which  gives  T  =  0.  The  cl  values  for  these  various  curves  are 

listed  in  Table  1.  The  final  value  of  a  considered,  a^.is  defined  so  that 

a,  -  a  is  5  percent  of  a.  -  d _ 

f  eq  r  l  eq 

Figure  1  demonstrates  several  interesting  features.  First,  the 
line  of  Equation  (10)  generally  falls  between  the  dissociation  and  recom - 
bination  curves  for  any  given  set  of  conditions  It  eliminates  the  variation 
caused  by  the  changing  parameters  and  leaves  only  the  variation  between 
the  dissociation  and  recombination  curves.  Although  the  factor  t  does  not 
bring  all  the  curves  together  it  would  appear  to  be  as  useful  a  single 
parameter  as  can  be  expected  for  correlating  all  these  curves.  Actually, 
knowledge  of  the  value  of  t  allows  one  to  estimate  the  reaction  time  only 
within  a  factor  of  10.  For  the  larger  values  of  AR0/BP  the  curves  depart 
from  the  logarithmic  part  much  sooner  but  the  recombination  rate  follows 
the  logarithmic  curve  closer  than  the  dissociation  curve. 

More  accurate  values  of  the  time  for  any  given  reaction  may  be  found 
by  using  the  appropriate  value  of  AR0/Bp  to  give  the  time  between  any  two 
values  of  a.  The  correct  value  of  a^  for  which  the  curve  are  drawn  must  be 
us  ed  in  making  the  computation. 

In  Figure  1,  all  of  the  curves  show  a  higher  dissociation  rate  than 
recombination  except  the  one  for  h/C^B  =  0.1  and  AR0/Bp  =  10^.  The 
suggested  conclus ion  is  that  for  the  value  of  parameters  of  most  practical 
interest  for  oxygen  and  nitrogen,  dissociation  is  always  faster  but  that  the 
opposite  condition  may  occur  at  values  of  AR9/Bp  lower  than  10^,  The 
relative  speed  of  the  processes  will  be  considered  in  more  detail  later. 

The  range  of  variation  of  the  parameter  t  is  shown  in  Figure  2  for 

the  same  range  as  the  parameters  considered  in  Figure  1.  r  varies 
7  14 

between  10  and  10  over  this  range.  The  spread  of  a  factor  10  in  ( t^  - 1) /r 
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shown  in  Figure  1  is  only  a  small  part  of  the  total  change  in  reaction  rate 
and  most  of  the  change  is  accounted  for  by  the  characteristic  time  t..  The 
value  of  A^/B0w  is  about  3.2  (10*^)  sec'  *  for  and  1.4  (10*4)  sec"^  for 
Oz  so  r  varies  from  the  order  of  1  sec  for  AR0/BP  =  10*®  to  the  order 
of  10'7  for  AR0/BP  =  106. 

An  examination  of  the  relative  effective  rates  of  dissociation  and 
recombination  can  best  be  made  by  considering  the  second  derivative  of 
Equation  (7).  A  positive  value  of  this  derivative  would  mean  that  the 
dissociation  rate  would  become  larger  than  the  r ec ombination  rate  for 
small  departures  from  equilibrium  and  a  negative  value  would  mean  the 
opposite.  If  this  derivative  is  taken  and  evaluated  at  equilibrium  it  is 


A  calculation  of  this  derivative  over  the  range  of  parameters  being 
considered,  gives  the  results  shown  in  Figure  3.  This  figure  shows  that 
the  derivative  is  positive  over  most  of  this  range  and  values  down  to 
AR0/Bp  =  104  must  be  considered  in  order  that  the  derivative  is  negative 
over  a  large  range  of  h/C  0.  These  lower  values  of  AR0/Bp  correspond 

v4 

to  high  pressures,  about  10  atm  for  N^.  While  there  seems  to  be  no 
fundamental  reason  why  dissociation  should  be  faster  than  recombination, 
the  parameters  for  hypersonic  flight  in  air  are  such  that  they  fall  in  the 
zone  where  this  is  the  case. 
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3.  DIRECTION  OF  REACTIONS  IN  ISENTROPIC  FLOW 

In  general  in  an  adiabatic  flow  field  recombination  will  occur  through 
an  isentropic  expansion  of  the  gas.  However,  it  should  be  recognized  that 
the  reaction  is  not  always  in  this  direction,  while  a  decrease  in  temperature 
tends  to  decrease  the  dissociation,  a  decrease  in  density  tends  to  increase 
it.  The  expression  for  a  change  in  entropy  as  a  function  of  T  and  a  is  given 
by  Lighthill 


dS  = 


C  dT 
v 

RT 


(1 


(17) 


If  the  equilibrium  expression  for  a  as  a  function  of  p  and  T, 
Equation(8),  is  used  to  eliminate  p  and  dS  is  set  equal  to  zero,  then 


C 

(1  +  u)  0/T  -  -j£- 

e/r  +  (i  +  a)(|  * 


The  condition  for 


da 

dT 


0  is  (1  + 


a) 


(18) 


(19) 


For  Cv/R  =  3  a  plot  Of  this  condition  is  shown  in  Figure  4  on  the  T,  p  plane. 
For  all  conditions  above  and  to  the  right  of  this  curve  dissociation  will 
occur  through  an  equilibrium  isentropic  expansion  and  for  cases  below  and 
to  the  left  recombination  will  occur.  Almost  all  cases  of  practical  interest 
lie  below  and  to  the  left  so  that  recombination  should  occur  through  an 
isentropic  expansion. 


4.  EFFECT  OF  CHEMICAL  NON- EQUILIBRIUM  OF  THE  FLOW  FIELD 

The  conditions  in  the  flow  field  cause  the  chemical  reactions  to  take 
place  as  has  been  illustrated  and  the  chemical  reactions  have  an  effect 
upon  the  flow  field.  A  chemical  reaction  requires  a  driving  potential  which 
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generally  leads  to  an  irreversible  change.  In  the  two  limits  of  equilibrium, 
no  driving  potential,  and  frozen,  no  reaction,  the  process  is  truely  revers¬ 
ible.  Since  an  actual  case  can  only  approximate  these  limits,  it  must 
always  be  somewhat  irreversible,  While  a  frozen  process  is  reversible 
while  it  reminas  frozen,  in  certain  ways  this  is  an  unrealistic  case.  For 
instance,  an  expansion  about  a  body  may  be  frozen  and  reversible  but  in 
the  wake  of  the  body  the  relaxation  will  eventually  take  place  in  an  irre¬ 
versible  manner.  For  the  equilibrium  case  there  will  be  no  subsequent 
relaxation.  In  this  respect  the  reversible  frozen  expansion  prepared  the 
way  for  the  irreversible  relaxation  unless  there  is  some  intermediate  step 
such  as  a  frozen  compression  which  restores  equilibrium  conditions 
reversibly. 

In  the  case  of  flow  about  bodies  consisting  of  a  rapid  pressure  rise 
through  a  bow  shock  and  then  an  expansion  to  ambient,  three  different 
chemical  rates  give  the  limiting  conditions  on  the  wake  flow.  Rapid  chem¬ 
ical  rates  give  equilibrium  conditions  throughout,  intermediate  rates  give 
an  equilibrium  shock  and  a  frozen  expansion  and  relaxation  in  the  wake, 
and  very  slow  rates  give  frozen  conditions  throughout  and  no  need  for  an 
eventual  relaxation  if  the  temperature  behind  the  expansion  is  not  high 
enough  to  cause  dissociation.  The  ratio  between  the  free  stream  velocity 
and  the  wake  velocity  for  nitrogen  under  these  three  conditions  and  several 
free  stream  Mach  numbers  are  as  follows: 


1)  2)  3) 


M 

Equil.  Shock 
Equil.  Expansion 

Equil.  Shock 
Frozen  Expansion 

Frozen  Shock 

Frozen  Expansion 

i 

15 

0.  899 

0.  765 

0.  904 

k 

20 

0.  794 

0.  656 

0.  918 

25 

0.  781 

0.  583 

0.  929 
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The  additional  decrease  in  the  wake  velocity  caused  by  the  frozen 
expansion  is  of  the  same  order  of  magnitude  as  the  decrease  in  velocity 
caused  by  the  normal  shock  for  the  equilibrium  case.  The  irreversibilities 
associated  with  the  relaxation  following  the  frozen  expansion  are  of  about 
the  same  order  as  those  introduced  by  the  equilibrium  normal  shock.  In 
other  words,  the  irreversible  chemistry  can  be  as  important  in  producing  the 
inviscid  wake  as  the  normal  shock.  The  case  of  frozen  shock  frozen 
expansion  gives  the  highest  wake  velocities  because  the  irreversibility 
associated  with  the  frozen  shock  wave  is  less  than  for  the  equilibrium 
shock  since  the  modes  of  energy  absorption  are  reduced  and  higher  tem¬ 
peratures  are  realized.  Since  the  frozen  shock  is  followed  by  a  frozen 
expansion  there  is  no  irreversible  relaxation  process. 

The  fact  that  the  ratio  of  the  wake  velocity  to  the  free  stream 
velocity  goes  up  with  Mach  number  for  the  frozen  shock,  at  first  seems 
surprising.  In  fact  this  velocity  ratio  approaches  l  as  Mach  number 
approaches  oo  if  the  specific  heat  of  the  gas  remains  finite.  A  velocity 
ratio  of  one  does  not  mean  that  there  are  no  total  head  losses  or  that  the 
temperature  in  the  wake  is  the  same  as  the  free  stream  temperature. 

As  M  — >  oo,  the  velocities  relative  to  the  object  has  become  so  large  that 
the  change  in  the  wake  velocity,  although  it  is  a  large  number  compared 
with  free  stream  speed  of  sound,  is  negligible  compared  with  the  velocity 
relative  to  the  object. 


5.  CONCLUSIONS 

a.  A  characteristic  time  can  be  derived  which  specified  the  reaction 
rate  near  equilibrium  and  gives  the  correct  order  of  magnitude  even  far 
from  equilibrium. 

b.  The  speed  of  the  dissociation  and  recombination  reactions  are 
the  same  near  equilibrium  and  their  relative  magnitudes  far  from  equilib¬ 
rium  depends  on  the  parameters  of  the  reaction.  For  the  normal  aero¬ 
dynamic  range  the  dissociation  is  the  faster, 

c.  Nonequilibrium  chemical  relaxation  can  introduce  irreversi¬ 
bilities  of  the  same  order  as  a  normal  shock. 
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Table  1.  Conditions  for  Reaction  Rate  Curves  Shown  in  Figure  It 


AR0 

W 

h 

Cv0 

d 

eq 

Dissociation 

a.  d, 

l  f 

Recombination 

a .  a, 

i  f 

iO6 

0.  1 

0.  04502 

0 

0.  04277 

0.  300 

0.  05777 

0.  3 

0.  49628 

0 

0.  4714 

0.  900 

0. 51647 

0.  5 

0.  92309 

0 

0.  87694 

1,  0 

0,  9269 

10® 

0.  1 

0,  08510 

0 

0. 50850 

0.  300 

0.  9585 

0.  3 

0.  59266 

0 

0.  5630 

0.  900 

0.  6080 

0.  5 

0.  99602 

0 

0.  94621 

1.  0 

0.  9962 

101® 

0.  1 

0.  11900 

0 

0,  11305 

0.  300 

0.  1280 

0*  2 

0.  38479 

0 

0.  3655 

0,  600 

0.  3956 

0.4 

0.  91386 

0 

0.  8681 

1,  0 

0.  9181 
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1* 

a,b, 

G 

2. 


3, 


4. 


Table 


Comparison  of  actual  reaction  tiroes  with  those  predicted  by 
logarithmic  law. 

Characteristic  reaction  time  as  a  function  of  reaction  conditions 
parameters. 

Second  derivative  of  reaction  rate  at  equilibrium  as  a  function  of 
reaction  conditions  parameters. 

Areas  of  reaction  conditions  for  which  dissociation  and  recombination 
occur  for  equilibrium  isentropic  expansion  of  a  diatomic  gas. 


i.  Conditions  for  reaction  rate  curves  shown  in  Figure  i. 


Figure  la.  Comparison  Of  Actual  Reaction  Times  With 
Those  Predicted  Ry  logarithmic  law 


Figure  lc. 


Figure  3.  Second  Derivative  Of  Reaction  Rate  At  Equilibrium 
As  A  function  Of  Reaction  Conditions  fturameters 
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